Reconstruction of f(R), f(T) and f(Q) models inspired by variable 

deceleration parameter. 
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ABSTRACT 

We study an special law for the deceleration parameter, recently proposed 
by Akarsu and Dereli, in the context of f(R), f(T) and f{Q) theories of modi- 
fied gravity. This law covers the law of Berman for obtaining exact cosmological 
models to account for the current acceleration of the universe, and also gives the 
opportunity to generalize many of the dark energy models having better consis- 
tency with the cosmological observations. Our aim is to reconstruct the f(R), 
f(T) and f(Q) models inspired by this law of variable deceleration parameter. 
Such models may then exhibit better consistency with the cosmological observa- 
tions. 

Subject headings: Accelerating universe, modified f(R), f(T) and f(Q) gravities, 
deceleration parameter 
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Introduction 



It is strongly believed that our universe is now experiencing an accelerat ed expansion 



This belief is supportec 



1998 



2004 



Perlmutter et al. 



by the recent observatio ns from type la supernovae fR iess et al 



1999 



Astier et al. 



20061 ). Large Scale Structure (|Abazaiian et al 



20031 ). Dark energy 



20051 ). and Cosmic Microwave Background anisotropics (jSpergel 
was the first idea to explain this accelerated expansion in the context of general relativity. 
This energy is a concept that we usually use for the unknown energy source in general 
relativity that is believed to be responsible for the observed acceleration of the universe. 
Modified gravities, on the other hand, are alternative wa ys for explaining this accele ration. 



Among the common f(R) modified theories of gravity (INoiiri and Odintsov 



Capozziello et al. 



2009 



Sotiriou and Faraoni 



2010 



Cruz and Dobado 



2007 



(Noiiri et al. 


2005 




2006 




Carter and IN 


eupane I 


2006 


Elizalde et al. 


201ol: 


Mvrzakulov et al. 


2011 


) which is 



20061 ) explaining 
propo sed 



Carter and Neupane II 



2006 



string effective action. In this proposal, the current acceleration of the universe is caused by a 
mixture of scalar phantom and (or) potential/stringy effects. On the other hand, a theory of 
f(T) gravity has recently been received attention. Models based on this modified teleparalle l 



gravi ty were presented 



20081 ) and dark ener gy models (|Bengochea and Ferraro 



in f(T) cosmology ( iJamil et al 



ternative to inflationary mod els (IFerraro and Fiorini 



20071 . 



2009). Moreover, attractor solutions 



20121 ). generalized second law of t hermodynamics in 



f(T) cosmology with power-law and 
power-law solutions in f(T) gravity ( 



logarithmic corrected ent ropies 



modified f(T) gravity ( IHoundjo et al. 



Setare and Darabi 



Bamba et al. 



20121 ). 



20121 ). and cylindrical solutions in 



20121 ) are among the interesting recent works in this 



context. A variant of models for modified gravities has also been obtained 



Noether symmetry 



2012 



Jamil et al. 



2011 



approach (ICapozziello and De Felice 



Wei et al. 



2012 



2008 



Vakili 



Atazadeh and Darabi 



20121 ). 



20081 : 



jy using the 



Hussain et al. 
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Recently, Akarsu and Dereli proposed a special law for the deceleration par ameter 



which is linear in time with a neg ative slope. This law covers the law of Berman (IBerman 



1983 



Berman and Gomide 



19881 ) (where the deceleration parameter is constant) used for 



obtaining exact cosmological models, in the context of dark energy, to account for the 



current acceleration of the universe. More recen tly, a comparison o 



standard ACDM cosmology has been presented (lAkarsu and Dereli 



this law with the 



20121 ). and also LRS 



Bianchi type- 1 cosmo 



studied (lAdhav 



ogical model with linearly varying deceleration parameter has been 
20111 ). According to this law, only the spatially closed and flat universes 
with cosmological fluid are allowed and the universe ends with a big-rip. In principle, this 
new law gives the opportunity to generalize many of these dark energy models having better 
consistency with the cosmological observations. 



The linearly varying 



( lAkarsu and Dereli 



20121 ) 



deceleration parameter of Akarsu and Dereli is defined by 



aa 
a z 



H 

IP 



-kt +771. — 1, 



where a(t) and H(t) are time dependent scale factor and Hubble parameter, respectively, k 
and m are positive constants and an over dot denotes the time derivative. Solving ([1]) for 
the scale factor, we obtain 



a{t) = axe^ 8 ^ 1111 ^- 1 ), k > 0, m > 1, 
a(t) = a 2 (mt + c 2 ) 1/m , k = 0, m > 0, 
a(t) = a 3 e C3i , k = 0, m = 0, 



(2) 
(3) 
(4) 



where ai,a 2 ,a^ and 02,03 are constants of integration. In this paper, we aim to find the 
f(R), f(T) and f{G) models inspired by the above law of variable deceleration parameter. 
This is called reconstruction method. This method has already been used by people in 
different models. In fact, considering the modified f(R),f(T) and f(Q) gravity models as 
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effective descriptions of the underlying theory of dark energy, it is interesting to study how 
these modified gravities can describe the dark energy. On the other hand, according to the 
law of variable deceleration parameter, one may take the opport unity to generalize many of 



20121 ). 



these dark energy models with better observational consistency (lAkarsu and Dereli 
Therefore, motivated by the above comments on the importance of modified gravities 
and the law of variable deceleration parameter, it is very appealing to study how these 
modified gravities can effectively describe the law of variable deceleration parameter. 
Reconstruction for f(R) was first developed by N ojiri et al for gravity consi s tent with 



reali stic cosmology and universe expansion history ( INojiri and Odintsov 



2005 



Nojiri et al. 



200 61). Another var iant of reconstruction is developed by Nojiri et al in terms of e-folding 



( INojiri et al. 



Recently, reconstruction of f(T) gravity is also studied in ( IBamba et al. 



singularity structure and it 


(Noiiri and Oc 


intsov 


2011) 


IBamba et al. 


2012 


)• 



In section 2, we study the above three cases in the context of modified f(R) gravity 
to find the corresponding f(R) models. In section 3, we study these cases in the context 
of teleparallel equivalent of General Relativity to find the corresponding f(T) models. In 
section 4, we study these cases in the context of scalar- Gauss- Bonnet gravity to find the 
corresponding f(G) models. The paper ends with a brief conclusion. 



2. Variable Deceleration Parameter and f(R) Model 



The most popular theory of modified gravities is the one kno wn as f(R) theory of gravity . 



'he action for this theory 



(ICapozziello et al. 



201 



2009 



coupled with matter S m is given by flNoiiri and Odintsov 



Capozziello and Francaviglia 



2007). 



20081 1 and ( iSotiriou and Faraoni 



S 



1 

2^2 



d 4 x^/^gf(R) + S m , 



(5) 
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where k 2 = 8nG. In agreement with the current observations, we take the spatially-flat 
Friedmann- Robert son- Walker (FRW) metric 



ds 2 = -dt 2 + a(t) 2 ^(dx 1 ) 2 . 



i=l 



Moreover, we assume the matter source to be a perfect fluid. By using the metric 
the perfect fluid in the Lagrangian (J5J) we obtain the field equations 

f(R) 



f(R) 



3(H 2 + H) f R (R) - 18 AH 2 H + HH f RR (R) + n 2 p = 



H + 3H 2 )f R (R) + 6 [8H 2 H + AH 2 + 6HH + dH/dt) /rr(R) 



(6) 
and 

(7) 



+ 36UHH + H) f R R R (R) + K 2 P = 



where f R = f'(R), f RR = f"(R), f RRR = f'"(R), the Hubble rate H is defined by H = a/a 
and the scalar curvature R is given by 



R = 12H 2 + QH. 



(9) 



Here, p and p are the energy density and pressure of the matter source, respectively. The 
equation of energy conservation is also obtained as usual 



p + 3H(p + p) =0. 



(10) 



Case 1 

By using Eq.(j2]) in Eq. flTUl) . we obtain 

p(t) = C x 



t 



3(l+ui) 
m 

kt — 2m 

where use has been made of the barotropic equation of state p 



wp with w being a 



constant, and C\ is a constant of integration. We then find 



Bit) = 2 km' 





'kt " 


2 


1 - 








m 





121 
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H(t) = 4k 2 m- 3 
Inverting H(t) to obtain the inverse function leads to 





~kt 1 


2 


-2 


"fct 


1 - 








1 




m 






m 



(13) 



Hm ± VH 2 m 2 -2 Eh 
Hk ' 



(14) 



where the positive sign will be selected so that t > 0. Moreover, by using Eq.fJTJ, we easily 
obtain 

H = -(q + l)H 2 . (15) 
Inserting this result in the right hand side of Eq.flH]) gives the following equation 

R = QH 2 (l-q). (16) 

By inserting t(H) in Eq. liTTj) . and using Eqs.Q, (fT3|) and (|T6"1) we get 

2 



3(l+u) 

p(R) = de^^k 



1 - 



l + ./l 



2fc 



m 2 ^fi/eCl-g) 



(17) 



The Friedmann equation (J7J) takes on the following form 



2R 2 (l + q) 
(l-q) 2 



V 6(1 -g) y y 6(l- 



q) 



f RR (R) 



f(R) 3qR 



6(1 -q) 



fn(R) 



2r< J-^ k 



+ k 2 de 



l+./l- 



2k 







m 2 y/R/6(l-q) . 

Note that g can be expressed in terms of R through Eqs. (jl4p . and (TT6l) . Therefore, in 
principle, the above equation is a differential equation for f(R). There is no known 
analytical solution for this equation, therefore one can look for numerical or approximate 
solutions. 
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Case 2 



Now, we use the second case for the scale factor. Putting Eq.Q in Eq.f JT0|) . and using 
H = (mt + c 2 ) _1 and Eq.(j9]) leads to a simpler form of p as 

p(t) = Arimt + c 2 )^-^ = ^V 3(1+w) = A 1 (H 2 ) 2m , 

«2 



or 



p(R) = M 



R 



3(1+™) 
2m 



(19) 



.6(1 -q) 

where A\ is a constant of integration and use has been made of Eq. (fl6l) . Inserting Eq.(fT9" 
in Eq.^7]) leads to the following differential equation for f(R) 

?2H , _\ / / 1 i ~\4/3 r>\ 3/2 



1T^F" 36 1 6(l- 9 ) J +18g 6TI^) 



/(it!) 3gi? 



i? 



3(l+w) 
2m 



0. 



(20) 



1.6(1 -g) J 

As in the previous case, there is no known analytical solution for this equation, so one can 
look for numerical or approximate solutions. 



Case 3 



For the third case, we have H = Const and R = Const. The Friedmann equation (J7J) 



reads 



f(R) 



+ 3H*f R (R) + k 2 Po = 



(21) 



where both H and p are constant. The solution of this equation is simply obtained as 

f{R)=A 2 exp(J^) +2k 2 p , (22) 



where A 2 is a constant of integration. 
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3. Variable Deceleration Parameter and f(T) Model 



The theory of modified gravity based on a modification of the teleparallel equivalent 
of General Relativity is calle d f(T) theory of gravity. The a ction for such a the ory coupled 



with matter L m is given by (IBengochea and Ferraro 



2011) 



2009|), flLinder 



20101 ) and flCai et al 



S = ^J d A xe[T + f(T)] + S m , 
where e = det{e 1 ^) = \J—g. The teleparallel Lagrangian T is defined as 



(23) 



where 



rp Q pprpp 

1 — °p 1 [IVl 



T% v = e£(0„e£ - d v e 
1 



(24) 



s p ^ = -(k^ p + 5^t%-5;t^), 



and is the contorsion tensor 



ZtT^V frppv rpup rp pu\ 

n p — 2^ p ~ p p I' 

The dynamical variable in the teleparallel equivalent of General Relativity is the vierbein 
e* . The field equations are then obtained 
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(25) 
(26) 

^ = ~2(l + / r - 12^/tt)' (2?) 
where fa = f'(T) and for = f"(T). Now, we investigate f(T) for three cases obtained in 
Eqs.©, (ED, and Q as follows. 



H 2 = ^ - \f - 2H*f T , 

6 

K 2 (p + P ) 
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Case 1 



In this case, Eqs.flTT}, 02}, UM and (E5J, lead to 



3(1+1^) 

p(T) = C 2 e m k 



1 + 



2fc 



(28) 



n 2 -v/-T/6_ 

where C*2 is a constant of integration and the Friedmann equation ( 126]) casts in the following 
differential equation for f(T) 



2Tf T - f(T) +T + 2K 2 C 2 e^^k 



l + ./l 



2k 



0. 



(29) 



Like some previous cases, there is no known analytical solution for this equation, and 
one can look for numerical or approximate solutions. 



Case 2 

Putting Eq.© in Eq.flEU, and using H = (mt + c 2 ) -1 and Eq.(j25H leads to 

P (t) = Mmt + c 2 )^ = = * (H") 8 ^ = * (-1) ^ , (30) 

where A 3 is a constant of integration. Inserting Eq. (|30|) in Eq. (j26|) leads to the following 
differential equation 

3(1+™) 
/ T*\ 2m 

2Tf T -f(T)+T + 2K 2 A 3 l--\ =0, (31) 
whose solution is obtained as 



4 3Q+w) „ 3(l+w) K 1 



l 2 2m 3 2m 3(1 + ™) 

f( T) = - T + B 1 VT + A 3 1 o3(1+m) 8 (-T)^, (32) 

1 2 2m 

where i?i is another constant of integration. 



Case 3 

Now, we use the third case for the scale factor. For this case, we have H = Const and 
T = Const. So, Eg. (1261) becomes as follows 



2Tf T -f(T) + T + 2K 2 p = 0, 



(33) 



where p = p = Const. The solution of this equation is simply obtained as 



f(T) = A 4 VT-T + 2K 2 p 0} 



(34) 



where A4 is a constant of integration. 



4. Variable Deceleration Parameter and f{Q) Model 



Now, we consider the following f(Q) action which descri bes Einstein's gravity 



with perfect fluid p lu s a function of the G auss-Bonnet term (INojiri and Qdintsov 



(INojiri et al. 



20061 ). (IRastkar et al. 



20121 ) 



coup 



20051 ). 



ed 



S= 2^J d ' X ^ ^ R + f^ + S " 



(35) 



where the Gauss-Bonnet invariant is defined by Q = R 2 - A R^R^ + R m x 9 W^ (J . The 



Friedmann equation for a flat FRW background is obtained (IRastkar et al. 

4# 2 + Qh - f(G) - 2AGH 3 f gg + p = 0, 
where p satisfies the conservation equation (fTQ|) . 



20121 ) 



(36) 



For the FRW metric, the Gauss-Bonnet term and Ricci scalar take the following forms 



g = 2A(HH 2 + H 4 



R = 6(H + 2H 2 ). 



(37) 
(38) 
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By using Eq.([T]), we easily obtain 

H = -(q + l)H 2 , (39) 

and inserting this in the right hand side of Q in Eq. (l37|) gives the following equations 

g = -2AqH\ (40) 

g = 2AH 4 [k + 4Hq{q + l)]. (41) 
Now, by using ( I4ip . we can write the Friedmann equation ( 136]) in the following form 

- - 2 H 2 + gf g - f(g) - (2A) 2 H 7 [k + AHq(q + l)]f gg + p = 0. (42) 

In this equation, p, q and H implicity are functions of g by following this algorithm: 
By using the function H(t), one may solve it for t = t(H). Then, it is possible to 
write H = H(H), and obtain H = H{g) by using Eq. (137j) . Thus, since p = p(H) and 
q = q(H,H), then we have p = p{g) and q = q(g). Therefore, we read the Eq.( ]4"2"]) as a 
differential equation for /(£7) which may be solved analytically or numerically. Now, we 
analyze the models given by (J2J), (j3J) and (T4]). 



Case 1 



By using Eqs.flTT}, fll2p, (TO, and MOj) we find 



p{g) = C 3 e~^k 



1+. 1- 



2k 

WW 



(43) 



where C3 is a constant of integration. By substituting the above p(G), and also H(g) from 
Eq. (T40l) . in Eq. (142l) we find the following differential equation for /((/) 

3 ( g n 1/2 



24g 



+gf g -f(g) 



(44) 



13 



(24) 



■2, Q\ 7,A 



24q) 



k + A 



24q 



1/4 



+ C 3 e m k 



1+. 1 



0. 



fgg 



2k 

WW 



No analytical solution for this equation is known, therefore one can look for numerical or 
approximate solutions. 



Case 2 



In this case we obtain 



whose time derivative becomes 



H 



mt + C2 



H = -mH 2 . 



Comparing the above result with the following equation derived from (1371) . namely 

H 



H\ 



2AH 2 

leads to the fourth order equation for H = H{Q) as follows 

G 



0. 



24(1 - m) 

Solving this fourth order equation is easy and gives the following solution 

H(G) = ± — -A r 

y J [24(1 -m) 

g ' 

24(1 -m) 

It is easy to obtain the following expression for p 



or 



H\g) 



1/4 



1/2 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 



p = p (a 2 cl /m r 3il+w) (l-mH) 



3(l + m) 



(51) 
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and by inserting p in Eq. (1421) . we obtain the following differential equation for f(Q) 

- 2 H 2 {Q) + Gf g -f{Q) (52) 
- {2AfH\g)[AH{G)m{m-l)]f gg 
+ p (a 2 cl /m )-^[l - mH{G)]^ = 0. 

Now, depending on the two choices for H{Q) in (149j) . the solutions of this differential 
equation can be obtained for /((?). 



Case 3 



For this case, we obtain H = Const and so a(t) = a 3 e . We may proceed as follows 

H = Const. =}► = H = -^—-H 2 , (53) 

which results in the following fourth order equation for H = H{Q) 

H*-^ = 0, (54) 
24 v ; 

with the solution 



H{G) = ± [y a ) ■ ( 55 ) 

For energy density we find p = po = Const, so by inserting H(Q) and p in Eq.(f42l) we obtain 
the following differential equation for f(Q) 

-^Uz) +Qfo-f(&)+Po = 0. (56) 
This equation is easily solved and we obtain the following solution 

f(g) = ^2-g^g + C& + Po , (57) 

where C4 is the constant of integration. 
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5. Final Remarks 

We have studied the recently proposed special law for the deceleration parameter by 
Akarsu and Dereli. According to this law only the spatially closed and flat universes are 
allowed and in both cases the cosmological fluid exhibits quintom like behavior, moreover 
the universe ends with a big-rip. This new law also gives the opportunity to generalize 
many of these dark energy models with better observational consistency. Motivated by this 
claim, we have tried to obtain the models of f(R), f(T) and f{Q) corresponding to this 
law in the framework of dark energy models. In this regard, we have obtained some exact 
solutions, and in some complicate cases we have left the solutions for numerical analysis. In 
a parallel way to that of Akarsu and Dereli, these solutions may give the opportunity to 
generalize many of the modified gravity models with better observational consistency. 
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